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We show how to calculate the magnetic-field and sheet-current distributions for a thin-film su- 
perconducting annular ring (inner radius a, outer radius b, and thickness d « a) when either 
the penetration depth obeys A < d/2 or, if A > d/2, the two-dimensional screening length obeys 
A = 2A 2 /d << a for the following cases: (a) magnetic flux $ z (a) trapped in the hole in the absence 
of an applied magnetic field, (b) zero magnetic flux in the hole when the ring is subjected to an 
applied magnetic field H a , and (c) focusing of magnetic flux into the hole when a magnetic field 
H a is applied but no net current flows around the ring. We use a similar method to calculate the 
magnetic-field and sheet-current distributions and magnetization loops for a thin, bulk-pinning-free 
superconducting disk (radius 6) containing a dome of magnetic flux of radius a when flux entry is 
impeded by a geometrical barrier. 



PACS numbers: 74.78.-w,74.25.Ha,74.25.0p,74.25.-q 
I. INTRODUCTION 

Recently Babaei Brojeny et ali reported exact ana- 
lytical solutions for the magnetic-field and sheet-current- 
density profiles for two current-carrying parallel copla- 
nar thin-film superconducting strips in a perpendicular 
magnetic field, fncluded were calculations for (a) the in- 
ductance per unit length when the two strips carry equal 
and opposite currents, (b) the zero- flux-quantum state 
when no net magnetic flux threads between the strips in 
a perpendicular applied field H a , and (c) the focusing of 
magnetic flux between the two strips in a field H a when 
each strip carries no net current. These problems are of 
relevance to the design of superconducting thin-film de- 
vices, especially superconducting quantum interference 
devices (SQUIDs). 

Of interest is the focusing of magnetic flux into the cen- 
tral hole in washer-type 2 SQUIDs and, in particular, the 
question of how much flux $h goes into the hole when the 
SQUID is in a perpendicular magnetic field H a = B a / fiQ 
and no net current circulates around the hole. The flux- 
focusing problem was examined by Ketchen et al.i who 
expressed $^ in terms of an effective pickup area of the 
hole, A e ff = §h/B a , which in general is larger than the 
actual area of the hole, Ah, but less than the area occu- 
pied by the washer, A w . Accounting only for azimuthal 
currents, they considered a washer of circular geometry 
(an annular ring) and derived a simple theoretical expres- 
sion for the effective area, A e ff » {S/TT 2 )Ah{A w /Ah) 1 ^ 2 , 
the theoretical approximations used being valid only for 
Ah « A w . Experiments on a series of square wash- 
ers with A w /Ah up to 10 4 yielded results in excel- 
lent qualitative agreement with the prediction, but with 



A efS » l.lA h (A w /A h )^ 2 . 

Experiments by Dantsker et ali on SQUIDs made with 
narrow superconducting lines separated by slots or holes 
(for trapping flux quanta during cooldown in the earth's 
magnetic field) have revealed that the presence of slots 
or holes increases the effective area over the value for a 
solid washer. This effect was confirmed experimentally 
by Jansman et al.i who were able to account for the 
increased effective area by treating the slotted washers 
as parallel circuits of pickup inductances. 

In this paper we introduce an approach suitable for 
extension to calculations of the magnetic-field and sheet- 
current-density distributions in superconducting thin- 
film strips, rings, and narrow lines. We consider the 
idealized case for which the penetration depth A obeys 
A < d/2 or, if A > d/2, the two-dimensional screening 
length A = 2\ 2 /d obeys A << a, such that the key 
boundary condition is that the normal component of the 
magnetic induction is zero on the surface of the super- 
conductor. A complicating consequence is that the sheet- 
current distribution in the superconductor has inverse 
square-root singularities at the edges. While a mutual- 
inductance approach such as that used by Gilchrist and 
BrandtS, and Jansman et ali is always applicable, we 
show here that an approach taking into account the 
inverse-square-root singularities from the beginning is 
simpler and more efficient. 

The authors of Ref. 3, obtained the flux-focusing result 
by superposition. They first calculated the induced cur- 
rent flowing in the clockwise direction in an applied mag- 
netic induction B a assuming zero magnetic flux in the 
hole. They approximated this current using the known 
result for a superconducting disk with no central hole. 
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They next calculated the induced current flowing in the 
counterclockwise direction in the absence of an applied 
field assuming a given amount of magnetic flux <frjj in the 
hole. They approximated this current using the known 
result for an infinite superconducting sheet with a round 
hole in it. Finally they obtained the relation between 
B a and $/j by equating the magnitudes of the two cir- 
culating currents. In the present paper, we show how 
to calculate all properties without making the small-hole 
approximations used in Ref. 0. We show how to solve 
the flux-focusing problem directly, as well as by superpo- 
sition. 

Another problem of interest is the calculation of the 
magnetic-field and current-density distribution for the 
case of a bulk-pinning-free type-II superconducting disk 
of radius b and thickness d << b in which the entry of 
magnetic flux is impeded by a geometrical barrier J£ An- 
alytic solutions for the field and current distributions and 
the magnetization in strips subject to a geometrical bar- 
rier have been studied for the bulk-pinning-free case in 
Refs. and ^] and for the case of Bean-model bulk 
pinning (J c = const) in Ref. [n[ Numerical results for 
the field and current distributions and the magnetization 
in disks subject to both a geometrical barrier and bulk 
pinning with a -B-dependent J c have been presented in 
Ref. El In the following, we present an efficient method 
for calculating the field and current distributions and the 
magnetization in bulk-pinning-free disks subject to a ge- 
ometrical barrier. 

Our paper is organized as follows. In Sec. II, we out- 
line our approach and set down the basic equations. In 
Sec. Ill, we apply this approach to calculate the induc- 
tance of an annular ring of arbitrary inner radius. In Sec. 
IV, we calculate the current circulating around a ring re- 
maining in the zero-flux-quantum state while subjected 
to a perpendicular magnetic field. In Sec. V, we con- 
sider the flux-focusing problem and calculate the mag- 
netic flux contained in the center of a ring in an applied 
magnetic field when there is no net current around the 
ring. In Sec. VI, we calculate the magnetization loop 
for a bulk-pinning-free thin-film type-II superconducting 
disk subject to a geometrical barrier. We briefly discuss 
our results in Sec. VII. 



II. BASIC EQUATIONS 

We consider a thin-film superconducting annular ring 
in the plane z = 0, centered on the z axis, with inner and 
outer radii a and b and thickness d « a. We assume 
that either A < d/2 or A << a if A > d/2, as discussed 
in the introduction. By the Biot-Savart law, the z com- 
ponent of the magnetic field in the plane z = i a 1 ^^ 



H Z ( P ) = H a + ±-J^ G(j>,fJ)KM)dff, 



where H a is the applied field, -K^(p) is the sheet-current 
density in the counterclockwise direction, 

G(p, p') = K(k)/(p + p') E(k)/(p p'), (2) 



k = 2(pp') 1/2 /(p + p'), (3) 

and K and E are complete elliptic integrals of the first 
and second kind with modulus k. An important bound- 
ary condition we will use in this paper is that H z (p) = 
for a < p < b. The total current in the counterclockwise 
direction is 

/= fK^dp, (4) 

J a 

and the magnetic moment along the z direction is 



= 7T / p 2 K cj> (p)dp. 



(5) 



Another quantity of interest is the magnetic flux up 
through a circle of radius p in the plane z = 0ji&±4 



$ z {p) = poHairp 2 + G A (p,p')K^p'W, 



(6) 



where 



Ga( P , pf) = (P + P'W - k 2 )K(k) - 2E(k)] (7) 

and k is given in Eq. ©. 

In the following sections we present solutions of the 
above equations and determine the corresponding sheet- 
current density K${p) for four cases: (a) self-inductance 
L = Q z (a)/I when H a = 0, (b) the zero-flux-quantum 
state [$ z (a) — 0] in an applied field H a , (c) flux focusing 
in an applied field [calculation of <I> Z (a) when 1 = 0], and 
(d) geometrical-barrier effects in a thin disk of radius 
b containing a Lorentz-force-free magnetic-flux dome of 
radius a. In each case, we assume a spatial dependence 
of the reduced sheet-current density of the form 



4g(u) 



tiu^J (u 2 — a 2 ) (1 — u 2 ) 



(8) 



where u = p/b and a = a/b and g{u) is a polynomial 
containing N terms, 



JV 



g( u ) = Y] 9 m (- r) 

* * I — n 



rn — 1 



(9) 



(i) 



Although we are not certain that such a choice gives an 
exact solution in general, it reduces to known exact so- 
lutions in various limits [a — > 0, b — * oo, or (6— a) << b], 
all of which have inverse-square-root singularities at the 
sample edges. To determine the N coefficients, N — 1 
equations are obtained by setting H z (p n ) = 0, where 
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p n = a + n(b - a)/N and n = 1,2, ...,N - 1. The N- 
th equation depends on the case under consideration; for 
case (a) we use Eq. (@J for given /, for case (b) we use 
Eq. © and set <J?z(a) = 0, for case (c) we use Eq. |@} 
and set 1 = 0, and for case (d) we use Eq. © and set 
5 (5) = 0. 

For numerical evaluation of the integrals in Eqs. JQ, 
lH} , (J5J , and © , it is convenient to change variables using 

the substitution v = p'/b = \J a 2 + (1 — a 2 ) sin 2 </> and to 
define the functions 

m«) = 4 r /2 G ( w > wx^r" 1 " -2 ^, (io) 
in, = - / ! — ^r -1 w~ 2 #, ii 

7T Jo 1 - a 

4 W 2 ?; — 5 
/m = 1 / (^4)™-^, (12) 

M«) = - G A (u, v ){^—4)™-\- 2 d^ (13) 
7T J 1 - a 

and 

a„ m = h m (u n ), (14) 

where u n = p n /b = a + n(l — a)/N, and n = 1, 2, iV — 
1. For a < p < b (a < it < 1), Eqs. (JTUJl and (T5|) 
are principal- value integrals, evaluated by splitting the 4> 
integral into two parts, one from to $(u — e) and the 
other from $(it + e) to vr/2, where 

*(u) = sin- 1 y^Zj (15) 

and e is an infinitesimal. For the results presented here 
we have used e = 1CP 7 . 



III. INDUCTANCE OF AN ANNULAR RING 

To calculate the inductance, we set H a = in Eq. 
JTJ and define Ki$ = (Ii/b)Ki$, where the subscript I 
henceforth labels all quantities that are specific to cal- 
culations of the inductance. To evaluate the coefficients 
girn in 

N 

gi (u)=Y,9i m Cf^) m -\ (16) 
^-^ 1 — a 

m—l 

we use the TV equations 

AT 

OLInmgim = Pin, (17) 

m—l 
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FIG. 1: Reduced magnetic field Hi z = bHi z /Ii, reduced 
sheet-current density Ki^, and polynomial gi (multiplied by 
3) vs u = p/b obtained while calculating the inductance for a 
superconducting ring with a = a/b — 0.5. 



n = 1,2, ...,7V, where a Inm = a nm and /3/„ = for 
n < N, and aiNm = 4 and /3/jv = 1 for w = TV. These 
equations are obtained from Eqs. (JJJ, ijHJ, EJ, II10|) . and 
l|T3]l and H z (p n ) = for n < N, and from Eqs. (@J, JBJ), 
@, and (ITTJ for n = N. 

Numerical results for H/ 2 = WJj z ///, if/0, and gi vs 
it = p/6 for a = 6/2 (a = 0.5) are shown in Fig. 1. In this 
calculation, as well as in all others in Sees. Ill- VI, the 
magnitude of the reduced magnetic field for a < p < b 
was less than 10 -5 except for p very close to a or b, where 
the numerical results for the principal-value integrals in 
Eq. l|lUfl became less accurate. Results for gi m vs 5 are 
shown in Fig. 2. The inductance is calculated from 

L = * Is /I 1 =i* ) b$ lz (a), (18) 

where 

N 

= < Pm(u)gi m , (19) 
m—l 

and is shown in Fig. 3 as a function of 5 = a/b. Dashed 
lines in Fig. 3 show expressions valid in the limits of small 
and large 5: For 5 << 1, the inductance approaches 
Lq = 2poa [or $i z (a) — 25], as obtained by Ketchen at 
al., 3 and for 5 — ► 1, the inductance approaches 

Li = p R[M 8R / w ) - ( 2 - ln4 )]; ( 20 ) 

as obtained by Brandt 15 for a superconducting annulus 
of mean radius R and width w « R. [Here R — (a + 
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FIG. 2: Coefficients gi, n in the poiynomiai of Eq. 1161 vs 
a — a/b obtained whiie calculating the inductance of a super- 
conducting ring. 



FIG. 3: Reduced inductance of a superconducting ring, L/pob 
vs a — a/b, calculated from Eqs. 11811 and 119H . Dashed curves 
show approximations valid in the limits a — > and a — * 1. 



6)/2 = 6(1 + a)/2 and to = (6 - a) = 6(1 - 5).] Equation 
(|20[1 can be obtained fromi£. L = fi R[ln{8R/r) - 2], the 
inductance of a superconducting ring of radius R and wire 
radius r << R, by replacing r by w/4H The empirical 
formula 

L 2 = /i O 6[a-0.197a 2 -0.0315 6 + (l + a)tanh" 1 a], (21) 

where a = a/b, fits our numerical results for L within 
0.06%, and a plot of it is indistinguishable from the solid 
curve in Fig. 3. 

The magnetic moment associated with the circulating 
current can be calculated from 

N 

miz = Iinb 2 2J fm9im- (22) 

m— 1 



IV. ZERO-FLUX-QUANTUM STATE 

Consider an annular ring that has been cooled into the 
superconducting state in the absence of a magnetic field, 
such that no magnetic flux is trapped anywhere in the 
ring. When a perpendicular magnetic field H a is applied, 
a circulating current is induced, but the ring remains in 
the Meissner state, and the magnetic flux up through the 
hole remains zero (there are no flux quanta in the hole). 
The induced sheet current density is Kz<p = H a Kz4>, 
where the subscript Z henceforth labels all quantititcs 
that are specific to calculations for the zero-flux-quantum 
state. To evaluate the coefficients gzm in 

N 

9z{u)=Y.9ZmC^) m -\ (23) 



we use the N equations 

N 

^ OiZnmgZm = flzn, (24) 
m— 1 

n = 1,2,..., N, where a Znm = a nrn and (3 Zn = -1 for 
n < N, and azNm = 4>m{a) and /3zn = —na 2 for n = N. 
These equations are obtained from Eqs. QJ, (jHJ, JJjJ, IjlOfl . 
and (CH|) and H z (p n ) = for n < N, and from Eqs. I©, 
©, ©, and © for n = N. 

Numerical results for Hz z = Hzz/ H a , Kz<p, and gz 
vs u = p/b for a = 6/2 (a = 0.5) are shown in Fig. 4. 
Results for gzm vs a are shown in Fig. 5. The magnitude 
\Iz\ of the induced current is obtained from Iz = H a blz, 
where 

N 

Iz = ^ i m gzm, (25) 

m— 1 

and is shown in Fig. 6 as a function of a = a/b. Dashed 
lines in Fig. 6 show expressions valid in the limits of small 
and large a: For a << 1, the induced current approaches 
Iz = —^Hob/ft [or Iz = — 4/7r], as obtained by Ketchen 
at al.j* and for a — > 1, the induced current approaches 
Iz = -irR 2 B a /L 1 . 

V. FLUX FOCUSING 

We now solve for the current and field distribution 
when a superconducting annular ring is placed in a per- 
pendicular magnetic field H a subject to the condition 
that there is no net current circulating around the ring. 
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FIG. 4: Reduced magnetic field Hz z = Hz z /H a , reduced 
sheet-current density Kz<t>, and polynomial gz (multiplied by 
3) vs u = p/fe for the zero-flux-quantum state with a — a/b = 
0.5. 



FIG. 6: Magnitude of the reduced current, —Iz/bH a , vs a = 
a/b for the zero-flux-quantum state calculated from Eq. 1250 . 
Dashed curves show approximations valid in the limits a — > 
and 5 — ► 1. 
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FIG. 5: Coefficients gzm in the polynomial of Eq. I)2M|I vs 

a = a/b for the zero-flux-quantum state. 



We wish to determine how much magnetic flux is focused 
into the hole in the middle of the ring. The sheet current 
density in this case is Kp^ — H a Kp^, where the sub- 
script F henceforth labels all quantities that are specific 
to calculations of flux focusing. To evaluate the coeffi- 



cients gp m in 



JV 



9f 



1 — a 



we use the N equations 



N 

^ C*Fnm.gFr, 



(26) 



(27) 



n = 1,2,..., N, where a Fnm = a nm and Fn = -1 for 
n < N , and apNm — im and Pfn = for n = N. These 
equations are obtained from Eqs. (JIJ, JSJ, iJjJJ, lfTU|) . and 
lfl4")) and H z (p n ) = for rt < TV, and from Eqs. (@J, ©, 
@, mi, and / = for n = N. 

Numerical results for Hp z = Hp z / H a , Kp^, and gp 
vs u = p/b for a = 6/2 (a = 0.5) are shown in Fig. 7. 
Results for gp m vs a are shown in Fig. 8. The magnetic 
flux focused into the hole is <&p z {a) = p H a b 2 $p z (a), 
where 



$F Z (a) 



7T0 2 



E 



4>m{a)g Fr , 



(28) 



The effective area of the hole (which corresponds to 
the effective pickup area of a SQUID made of a circu- 
lar washer), defined via <&p z {a) = poH a A e f f , is always 
larger than the actual area of the hole, Ah — 7ra 2 . We 
find 



.4 



eff 



^Fz(a) 
p H a ira 2 



1 



N 



4>mia)gi 



(29) 



m— 1 
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FIG. 7: Reduced magnetic field Hfz = HFz/H a , reduced 
sheet-current density Kf<i>, and polynomial qf (multiplied by 
3) vs u = p/b for flux focusing with a — a/b = 0.5. 
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FIG. 9: Reduced effective area, A e ff/A h , vs 1/5 = fo/a for 
flux focusing calculated from Eq. I|29^ or (!■'-!( )t . Dashed curves 
show approximations valid in the limits 5-^0 and a — * 1. 



0.4 




0.0 0.2 0.4 0.6 0.8 1.0 
a/b 



FIG. 8: Coefficients gFm in the polynomial of Eq. (12(it vs 
a — a/b for flux focusing. 



which is shown in Fig. 9 as a function of 1/5, = b/a. 
Dashed lines in Fig. 9 show expressions valid in the lim- 
its of small and large a: For a « 1, A e ff/Ah ap- 
proaches (8/7r 2 )(6/a) [or 8/7T 2 a], as obtained by Ketchen 
at al.ji and for a — > 1, A e ff/A h approaches {R/a) 2 
[or((l + a)/2a) 2 ], where R — (a + b)/2 is the mean radius 
of the ring. 



The flux-focusing problem also can be solved from a 
linear superposition of the fields calculated in Sees. Ill 
and IV. From Kp<f, = Kj^ + Kz<f> and the condition Ip = 
I] + I Z = we obtain g F {u) = -I z gi(u) + g z {u), g Fm = 
—Izgim + gzm, and the result 

A h ira 2 
which gives numerically the same values as Eq. I|29|). 

VI. GEOMETRICAL BARRIER 

We next present an efficient method for calculating the 
magnetic-field and current-density distributions and the 
magnetization of a bulk-pinning-free type-II supercon- 
ducting disk subject to a geometrical barrier, which im- 
pedes the entry of vortices into the disk. We consider 
a disk (radius b and thickness d << b) in the plane 
z = 0, centered on the z axis, intially in the Meissner 
state. We assume that the London penetration depth 
obeys A < d/2 or, if A > d/2, that the two-dimensional 
screening length A = 2 A 2 /d obeys A << b. When a per- 
pendicular magnetic field H a is applied, a sheet-current 
densitjiii 

K,{p) = - A J±P= (31) 

7T yjb 2 - p 2 

is induced. The resulting magnetic field in the plane 
z = 0, determined from Eq. Q, is H z {p) — for p < b 
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andi& 



H z {p) = H a ll 



— = — sin ( - ) 



(32) 



for p > b. 

A geometrical barrier prevents vortices from entering 
the film until the magnetic field at the edge (accounting 
for demagnetizing effects) reaches the value H s . We ex- 
pect that H s = H c i, the lower critical field, if there is 
no Bean-Livingston barrier, or H s H c , the bulk ther- 
modynamic critical field, if the edge is without defects 
and thermal activation is negligible. An equivalent crite- 
rion is that the magnetic flux begins to penetrate when 
the magnitude of the sheet-current density at the edge 
reaches the value K s = 2H s m& To estimate H z or 
at the edge of the film, we note that the approximations 
that led to Eqs. (|31|l and (|32[l break down and that the 
inverse-square-root divergences in these equations are cut 
off when p is within 5 of the edge, where 8 is the larger of 
d/2 or A. Accordingly, we approximate H z at the edge of 
the film by replacing p in the square-root denominator of 
Eq. by b + 5 and using S « b, such that H z (edge) 
« (H a /ir)y / '2b/5. Similarly, we approximate at the 
edge of the film by replacing p in the square-root denom- 
inator of Eq. (|3f \ by b — 5 and using 5 « b, such that 
A^(edge) w — (2H a /n)y / 2b/S. Whichever criterion is 
used [H z (edge) — H s or jAT^edge)! = K s = 2H S ], we es- 
timate that the geometrical barrier is overcome when the 
applied field is equal to Hq = nH s y/8/2b. (In this paper 
we have chosen 5 — 5/b — 0.01, such that Ho = 0.222iJ s . 
See Fig. 13.) 

When H a > H such that H z (edge) > H s , vortices 
nucleate at the edge of the disk and move rapidly towards 
the center of the disk under the influence of the Lorentz 
force per unit length, f = J h x o , where Jh = VxH re „, 
(f>o is a vector of magnitude (f>o — h/2e along the vortex 
axis, and H rev is the thermodynamic magnetic field in 
equilibrium with the magnetic flux density B inside the 
superconductor. As more vortices enter, the return field 
outside the disk generated by the vortices inside the disk 
gradually reduces the value of the field at the edge to 
H s , thereby halting further vortex nucleation. If bulk 
pinning is ncglible, the case considered in this paper, the 
vortices adjust their positions such that the magnetic flux 
density (averaged over the intervortex distance) in the 
plane of the disk B z (p) has its maximum value at the 
center, decreases monotonically to zero at p = a, and 
remains zero for a < p < b. The corresponding sheet 
current density Kh4> — Jn^d is zero for p < a, such 
that the Lorentz force on any vortex vanishes and no 
further motion occurs. Screening supercurrents still flow, 
however, in the vortex-free region a < p < b. 

To good approximation when d << b, the result- 
ing magnetic-field and supercurrent distributions are the 
same as those generated by a thin superconducting an- 
nular ring (a < p < b) in a perpendicular applied field 
H a , when the solutions are subject to the constraint that 



the sheet current density is zero at p = a. The 
Biot-Savart law [Eq. and its extension to \z\ > 0] 
guarantees that the current density Jb = VxB/jUo is 
zero everywhere except within the ring a < p < 6; 
thus Kb<p = JB<pd is zero for p < a. Because Jh 
and Jb in thin films are dominated by the curvature 
of H re .„ and B//xo, rather than by the gradients \7H rev 
and VB//Jor^ it can be shown that the difference be- 
tween Kh<p and Kb<p is of order (d/b)H a , decreases for 
B > 2B c i as Hrev approaches B//io, and is negligible 
for the thin films considered in this paper (d/b = 0.01). 
Nevertheless, our simplified approach would be incapable 
of calculating details in the structure that has been ob- 
served in the magnetic flux-density distribution at the 
vortex-lattice melting transitional To treat such a prob- 
lem would require a more refined approach such as that 
in Refs. and I22I which calculates the local 3b cur- 
rents flowing at the vortex solid-liquid interface and dis- 
tinguishes between H rev and B//io. 

The magnetic-field and supercurrent distributions for 
the case of a thin pin-free disk subject to a geometrical 
barrier therefore can be calculated efficiently by using 
an approach similar to that used in Sees. II-V. When 
a Lorenz-force-free dome of magnetic flux occupies the 
region p < a, the sheet current density in the region a < 
p < b is Kg<p — H a KGti>, where the subscript G hencefore 
labels all quantities that are specific to the geometrical- 
barrier problem. To evaluate the coefficients gem in 



9g(u) 



N 

E 

m— 1 



,u, — u 

ga m {- — -) 



we use the N equations 



N 



OLFnm9Gr, 



PGn, 



(33) 



(34) 



n = 1,2,..., A, where a Gnm = a nm and (3 Gn = -1 for 
n < N, and acNm = Sim and /3gn = for n — N. These 
equations are obtained from Eqs. JHJl, JHJl, IjlOl) . and 
(|T1J| and H z (p n ) = for n < N, and from Eqs. ©, ©, 
and Kg4,{o) = for n = N. 

Numerical results for Hq z = Hq z / H ai Kg4>i an d 5g 
vs u — p/b for a — b/2 (a = 0.5) are shown in Fig. 
10. In these calculations, we have made no distinction 
between H rev and B//i , which corresponds to assuming 
that B w pqH. However, in cases for which B differs 
significantly from PqH, our plots of Hc z /H a (such as in 
Fig. HOfl would correspond most closely to plots of the 
reduced flux density BGz/B a . Results for gam vs a are 
shown in Fig. 11. The magnetic flux contained within 
p < a can be obtained from 



§Gz{a) = p H a b 2 $ Gz (d), 



(35) 



where 



N 

®Gz(a) = 7ra 2 + 2J <t>m(a)g Gm , (36) 

m— 1 
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FIG. 10: Reduced magnetic field Hqz = Hgz/Hc, reduced 
sheet- current density Kg^,, and polynomial go vs u — p/b for 
a pin-free disk of radius b with a geometrical barrier and a 
flux dome of reduced radius a = a/b = 0.5. 



FIG. 11: Coefficients gGm in the polynomial of Eq. 1331 vs a 
for a pin-free disk of radius b with a geometrical barrier and 
a flux dome of reduced radius a — a/b. We require gai = 0. 



and the average magnetic flux density in the disk is 
B av = <&Gz {a)/irb 2 . Figure 12 shows how B av /B a and 
Hgz(0)/ H a , where Hg z (0) is the magnetic field at the 
center of the disk, depend upon a. 

We next calculate the magnetization, i.e., the mag- 
netic moment per unit volume of the disk, Mq z = 
mGg/Ttvd, where mcz is calculated from Eq. (J5J. The 
initial magnetization of the disk in the Meissner state 
(0 < H a < H , see Fig. 13), calculated from Eq. ffity. 
\3m^ Mqz — ~XoH a , where xo — &b/3nd; i.e., the ex- 
ternal magnetic susceptibility^ in this case is \ — ~ Xo- 
Whenever there is a dome of magnetic flux within the 
region p < a, the magnetization, obtained from Eqs. JSJ, 
©, and may be calculated from 

3^ N 
M Gz = —XoH a ^2 fm9Gm, (37) 

m— 1 

where f m and gam depend implicitly upon a. 

For Ho < H a < Hi rr along the field-increasing mag- 
netization curve at the critical entry condition (see Fig. 
13), H a and a are related via 

N 

H a = -iWl-5 2 / 9Gm, (38) 

m— 1 

This equation follows from the condition that 
l-fTc^edge)! = 2H S , where Kc^iedge) is obtained 
by evaluating Eq. (JSJ at u = 1 but replacing y/1 — u 2 
in the denominator by \J25jb, as in the evalua- 
tion of Hq- When a — 0, we see by comparing 



Eqs. JSJl and lj3T|) that ga{ u ) = —u 3 , such that 
go™. = Sm4 (see also Fig. 11), / 4 = 8/37T, and 
MqzI — ~XoHo at H a = Ho- In the limit as a — ► 1, 
Kc<ji ~ — 2i/u — a/y/1 — u, such that gc( u ) ~ ~ ir{u — a), 
9Gm ~ -7r(l - a)5 m2 , /2 ~ 1, -ffa ~ H y/2/wy/l - a, and 
M Gz1 / Xo H s » -{3n 2 /8)5H s /H a , where 5 = £/b « 1. 

For i?o < ff a < -ffirr along the field-decreasing mag- 
netization curve at the critical exit condition, we assume 
that the radius a of the vortex-filled region has reached 
within S of the radius b of the disk; i.e., a = 1 — 6. Us- 
ing Eq. H37|) with ga m ~ —Tf8S m 2 and / 2 « 1, we obtain 
M GzL / X0 H s « -(3n 2 /8)5H a /H s . See Fig. 13. 

The field-increasing and field-decreasing magnetization 
curves in Fig. 13 meet at _ff a = i^ rr , the irreversibility 
field. The criteria we used for the critical entry and exit 
conditions lead to the result that Hi rr sa H s , where the 
magnetization is given by Mazirr/XoH s ~ — (37r 2 /8)i5. 
However, the above expressions for Hi rr , Mozim and 
MqzI are the least reliable results of our paper, because 
all these quantities are very sensitive to the precise con- 
ditions for entry and exit at the edge of the disk, in- 
cluding such details as the shape of the edgei2i2ii2ii2i2£ 
The magnetic moment responsible for the magnetization 
Mczirr and MqzI is produced by currents that flow only 
within a very narrow band around the disk's edge, where 
a theory more accurate than ours is needed. 

The minor hysteresis "loop," shown as the dashed 
curve in Fig. 13, can be calculated as follows. We start 
at a point on the field-increasing magnetization curve 
where the flux dome has radius a±. The magnetic flux 
contained within the dome <5>Gz(ai), the magnetization 
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0.0 0.5 1.0 

a/b 

FIG. 12: Reduced average flux density B av /B a (solid) and 
reduced flux density at the center Hcz{0)/H a (dashed) of a 
pin-free disk of radius b with a geometrical barrier and a flux 
dome of reduced radius a = a/b. 



H/H H/H H/H 

Q s l s iiT ; 




-0.25 I — lj — i — u — i — i — u — i — u — i — i — u — lj — i — u — i — i — I 
0.2 0.4 0.6 0.8 1 1.2 

H/H 

a s 

FIG. 13: Calculated hysteresis in the reduced magnetization 
Mgz/xoHs vs reduced applied field H a /H s for a pin-free disk 
with a geometrical barrier (solid). The dashed curve shows a 
reversible minor hysteresis "loop" occurring when the applied 
field is reduced after the applied field H a has reached Hi along 
the field-increasing magnetization curve. As H a decreases, the 
flux dome expands, but the flux contained within the dome 
remains constant. 



Mq z i, and the corresponding applied field Hi are ob- 
tained from Eqs. l|3"7|) . and where a, f m , and 
gem are all evaluated at a = di = ai/b. As the applied 
field H a is reduced from its starting value Hi , the radius 
a of the flux dome expands, but the magnetic flux within 
the dome remains constant. For each value of a > di, we 
recalculate f m , gam, and $Gz(a). We then use Eq. i|3"5|) 
to obtain the corresponding value of the applied field, 

H a = H 1 ® Gz {a 1 )/$ Gz {h), (39) 

and Eq. I|37|> to obtain the corresponding value of the 
magnetization. 

VII. DISCUSSION 

In this paper, we have presented an efficient method 
for the calculation of magnetic-field and current-density 
profiles for thin-film rings in the Meissner state and for 
bulk-pinning-free disks subject to a geometrical barrier. 
In each case, the sheet-current density was expressed in 
the form of Eq. JSJl, where the quantity g(u) in the nu- 
merator is a polynomical of degree N — 1. 

For all the calculations presented in the figures, for 
which we assumed N = 5, we found that the magnitude 
of <?5 was less than 0.0012 in each case (see Figs. 2, 5, 8, 
and 11) and that its contribution to g(u) was less than 
1.1%. Using N = 6 yields values of ge whose magnitudes 
are much smaller than those of 35, and the values of the 
calculated physical quantities are altered only in the sixth 
decimal place. 

Moreover, we offer the conjecture that the problems 
we solved numerically in Sees. Ill- VI might be solved 
analytically with functions gi,gz,9F, and go that are 
third-order polynomials in u\ i.e., the sums in Eqs. 0, 
(Unj, (E>J>> and might simply terminate with 

N = 4. As evidence in support of this conjecture, we 
note that our calculations for a — a/b — 0.1 and 0.5 
with N = 4, 5, 6, and 7 yielded values of L/^b [Eqs. 
(fT5|l and i(T§|) ] that differed only in the fifth decimal 
place. Similarly, values of Iz/H a b [Eq. A e ff/Ah 
[Eq. (J25J)], and M Gz] / X oH s [Eqs. © and Qg)] calcu- 
lated for d = a/b = 0.1 and 0.5 with N = 4, 5, 6, and 7 
differed at most only in the fourth significant figure. It is 
possible that the values we obtained for g^^ge, and gi in 
Sees. Ill- VI were nonzero only because of small numeri- 
cal errors introduced because we performed the integrals 
in Eqs. (|10|) - (|13|) numerically rather than analytically. 

Although in this paper we have considered only bulk- 
pinning-free thin-film rings and disks, it should be possi- 
ble to extend the present approach to develop an efficient 
method, complementary to that of Ref. 6, for numeri- 
cally calculating quasistatic magnetic-field and current- 
density distributions in rings and disks subject to both a 
geometrical barrier and bulk pinning. Such distributions 
recently have been calculated an alytically for infinitely 
long st rips in Refs. llll25l26l27l28L 
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